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Adiabatic quantum computation ( AQC) [U [2] , either in its universal form [3j S] , or in the form 
of adiabatic quantum optimization or quantum simulations [7|, presents a viable alternative 

to gate-model quantum computation (GMQC). Although a part of the original motivation for in- 
troduction of the AQC [5] was the promise of the increased stability against decoherence due to the 
energy gap between the ground and excited states, the question of the role of decoherence in AQC 
remains an open one. This uncertainty makes it important to quantify more precisely the decoher- 
ence properties of AQC. A crucial step towards this would be to define a quantitative characteristic 
of the decoherence strength in AQC, that plays a role similar to the decoherence time for GMQC. 
However, in the case of AQC, decoherence has qualitatively different, static effect on the qubits, not 
limiting the operation time of an algorithm [8]. In this work, we propose the ground state fidelity, 
defined as the distance between the open and closed system reduced density matrices normalized 
to the Boltzmann ground state probability, as a quantitative measure of decoherence-induced de- 
formation of the ground state in AQC, analogous to the decoherence time for GMQC. We calculate 
the fidelity perturbatively at finite temperatures and express it through the same environmental 
noise correlators that determine the decoherence times in GMQC. We discuss the relation between 
fidelity and the relaxation and dephasing times of the qubits, and its projected scaling properties 
with the number of qubits. 



In AQC, adiabatic evolution of the ground state of 
a qubit system realizes the solution of a computational 
problem represented by an appropriately designed Hamil- 
tonian, which is typically written as 

H s = A{s)H v + B{s)H P , (1) 

where s = t/tf with tf being the total evolution time. 
At s = 0, one has A(Q) = 1, B(Q) = 0, and the system 
is initialized in the ground state of the Hamiltonian Hjj, 
which usually consists of the uniform superposition of all 
computational basis states. The energy scales A(s) and 
B(s) are varied monotonically so that at s = 1, A(l) = 
and B(l) = 1. If the evolution is slow enough, an 
isolated qubit system stays in the ground state with high 
fidelity throughout the evolution, and at s = 1 reaches 
the ground state of Hp, which provides a solution to a 
computational problem. 

If the qubit system is weakly coupled to a dissipative 
environment, two effects are expected. First, the low- 
frequency part of the environmental noise moves the sys- 
tem energy levels relative to each other. This results 
in a dephasing of the energy eigenstates that eventually 
suppresses all off-diagonal elements of the qubit density 
matrix in the energy basis. However, since the popula- 
tion of the ground state is the only important part of the 
computation and the relative phases of the energy eigen- 
states do not carry any information, this does not affect 
AQC. Second effect of the coupling to the environment is 
that it induces the thermal transitions between the qubit 
energy levels pushing the qubit system towards thermal 
equilibrium at a temperature T. For slow evolution, the 
instantaneous probability to be in the ground state is 
given approximately by the Boltzmann distribution, and 



so the qubit system loses some of the ground state prob- 
ability due to thermal occupation of the excited states. 
Such a thermal loss of probability can be compensated 
by multiple iterations of an AQC algorithm as long as it 
does not scale exponentially with the size of the system, 
i.e. as log as the number of excited states within roughly 
the energy kpT above the ground state does not grow 
exponentially. 

The preceding arguments provide an intuitive expla- 
nation for the predicted robustness of AQC against local 
environmental noise in the limit of weak coupling [5Hl5j. 
When the strength of the coupling to the environment is 
increased without changing cither the Hamiltonian or the 
temperature, the qubit Boltzmann distribution is still not 
directly affected. However, it is known that the decoher- 
ence time of the qubits decreases with increased coupling, 
and strong coupling to the environment eventually makes 
the qubits completely incoherent, rendering them useless 
for quantum computation. In GMQC, qubit decoherence 
leads to computation errors which, without error correc- 
tion, completely destroy the computation process. This 
is why the qubits' quality factor, which is the ratio of the 
decoherence time and the gate operation time, provides a 
good measure of the qubit performance in GMQC. It is, 
however, unclear how an increase in coupling to the envi- 
ronment, or equivalently decrease in qubit quality factor, 
affects AQC. 

In this paper, we look closely at what happens to the 
eigenstates of the qubit system in AQC when coupling to 
the environment is non- negligible. To ensure consistent 
notation throughout this paper, symbols with (without) 
"~" denote quantities related to the coupled (uncoupled) 
qubit system and environment. We use letters m,n to 
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enumerate the eigenstates and eigenvalues of the qubits 
(e.g., \ri), E n ), letters v, fj, to enumerate the eigenstates 
and eigenvalues of the environmental degrees of freedom, 
and letters a, b to enumerate the eigenstates and eigen- 
values of the total system (qubits+environment). The 
total Hamiltonian is H = H$ + Hb + Hi, where Hb 
and Hi are the environment and interaction Hamilto- 
nians, respectively. In the absence of coupling, Hi=0, 
and the eigenstates of the total system are \a) = \n)®\v) 
with eigenvalues E a = E n + E v . When Hi ^ 0, the 
new eigenstates are \a), which typically are entangled 
superpositions of the unperturbed states \a). For weak 
coupling, | a) is very close to |a) and the effect of the envi- 
ronment is thermalization of the qubit system. Once the 
environment is averaged out, the equilibrium of the total 
system gives the Boltzmann distribution for the qubits: 



P„ 



E 



-(E n +E u )/T 
ZsZb 



-E„/T 



(2) 



where Z s = J2 n e ~ E " /T and z b = E,v e ~ Eu/T are the 
partition functions of the qubit system and the environ- 
ment. 

As the coupling increases, the deviation of \a) from \a) 
grows. In equilibrium, the density matrix of the total 
system still has the Boltzmann form psb = Y]„ Pa\o) (a|, 

where P a = e~ S ^ T /Z SB , with Z SB = E a e ' Sa/T bein S 
the partition function of the total system. However, the 
reduced density matrix ps — Ttb[psb] of the qubit sys- 
tem alone is no longer given by the Boltzmann distribu- 
tion. The deviation from the Boltzmann form provides 
a good qualitative measure of how strongly the eigen- 
states | a) are deformed in comparison to the unperturbed 
states. 

Quantitatively, we define the ground state fidelity as 
the Uhlmann fidelity [16] between the reduced density 
matrix ps and the "ideal" ground state density matrix 
Po = |0)(0|, normalized to the Boltzmann ground state 
probability Pq: 



F = Pn 



-1/2 



Tr 



PaPsyfpa = V Pa/ Pa 



(3) 



where Po = (0|/?s|0) is the equilibrium probability for the 
qubits to be in the ideal ground state when coupled to 
the environment. In the weak-coupling limit, no defor- 
mation of the eigenstates is expected. Then Po = Po, and 
Eq. @ gives F = 1. This shows that Eq. ^ correctly 
separates the effect of the quantum deformation of the 
ground state, which can be viewed as the result of vir- 
tual transitions to the excited states, from the thermal 
loss of probability. Qualitatively, the effect of the vir- 
tual transitions, expressed in P, is different from that of 
the thermal transitions in two important aspects. First, 
it persists even at T = 0, when all the thermal transi- 
tions are suppressed. Second, it depends on the strength 
of coupling to environment (or decoherence time of the 
qubits), while thermal equilibrium probabilities only de- 



pend on the energy eigenvalues and temperature. Never- 
theless, similarly to thermal transitions, the virtual tran- 
sitions reduce the occupation probability of the ground 
state by transferring it to the higher-energy states. 

We calculate the fidelity ^ perturbatively and relate 
it to measurable parameters of the qubit system and en- 
vironment. As appropriate for AQC, we assume that the 
coupling Hi is weak. This allows us to employ the per- 
turbation theory in Hi around the non-interacting state 
of the qubit system and the environment with both of 
them in equilibrium at the same temperature T, i.e., den- 
sity matrices p s = J2 n p n\ n )( n \ and Pb = Yl v p A u ){ v \-> 
where P n and P„ are the Boltzmann probabilities. The 
reduction of the ground state probability due to finite Hi 
is caused by a change 5P in the equilibrium probability 
Po as a result of renormalization of the energy eigenval- 
ues (Lamb shifts), and probability transfers into and out 
of the ground state due to renormalization of the wave- 
functions. The probability Po that defines the fidelity ^ 
can be expressed as 



P = Tr B ,s 



|0)<0|$> a |a)(a| 



(4) 



Introducing interaction-induced corrections to the equi- 
librium probabilities P a = P a + SP a , where P a — P n Pu, 
and wavef unctions: \a(n,v)) = \n)®\v) + \Sa(n,v)), we 
can rewrite this expression to the lowest non-vanishing 
order in Hi as 

P =SPa + J2 P n T ^B,s [|0)<0|®p B • \a(n, v)) (a(n, u)\] . 

n 

Using the relation |0)(0| = l-^ m J \m)(m\ to transform 
the n — term in Eq. ([5| we obtain 

p = Po + SPa - ]T (r „p - r n0 p„) , (6) 

where 

r m „ = (n\Tr B [\$a(m, is))(5a{m, v)\p B ]\n) . 

The terms proportional to T in Eq. ^ describe the reduc- 
tion of the ground state probability as a result of renor- 
malization of the qubit system wavefunctions by their 
interaction with the environment. 

Next, we calculate SPq and T mn . Quite generally, the 
interaction Hamiltonian Hi is 



(7) 



where a" are the Pauli matrices for the jth qubit, a = 
x,y,z, and q" are the corresponding operators of the 
noise generated by the environment. As usual, the aver- 
ages of the noise operators vanish, (g") = 0. Then, in the 
weak coupling regime, the effect of environment is fully 
characterized by the noise spectral densities: 



Sf(u)= / dte iut (qmqm) 



(8) 
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where (...) = TtbIpb---} is the average over the envi- 
ronmental degrees of freedom. For simplicity, we limit 
our discussion to the most typical case when the noises 
with different a or j are uncorrelated. It is shown in 
the supplementary information (SI) that the perturba- 
tion expansion in Hi in this situation gives 



8P = -pp Q ]T Gp„-MK«J 2 



r m „ — ^2Wj,nm\ 2 J 2 n 



du Sf(u) 



du Sf(u) 
2tt w mn + to ' 

(9) 



where <J° ! nm = (n\a"\m) and u) nm = E n —E m . Substitut- 



ing ([9| into Q and then into we obtain 

dw Sf(uj)(P n -8 n0 ) 

u j,nm\ I 



j,a,n>0 



duS?(u)-(P n /P )Sf(-w) 



4tt 



Equation ( 10 ) is our main result. It is well-defined at T = 
0, when all thermal excitations are suppressed, i.e., P n — 
for n > and = at cj < 0. Hence, the values of 

u) around — cj m o, when the denominator in (10) vanishes, 



do not contribute to the integral. When T ^ 0, the 
divergences that appear at oj = —oJ m o reflect the fact that 
environment can also create real thermal excitations of 
the qubit system. However, the detailed balance relation, 
= e - ^" Sf(uj), ensures that these divergences 



Sf(- 



cancel each other out and Eq. ( 10 ) is well-defined also at 
T ^ (see the SI). 



Equation ( 10 ) is now applied to specific problems. The 



first example we consider is a typical individual qubit with 
the Hamiltonian 



Ha 



-[ea z + Aa x }/2 



coupled as in Eq. ([7 1, but only through a z , to the envi- 
ronmental noise with spectral density S(ui) (J8|. In the 
usual weak-coupling approximation [see, e.g., Ref. [17], 
the qubit decoherence time T 2 * is given by 



T* 2Ti T v 



(12) 



where T\ and T v are the relaxation and pure dephasing 
times, given by 



T-f 1 = (A 2 /n 2 )[s(n) + s{-n)}, 

T- 1 = (e 2 /n 2 )S(0), 



(13) 
(14) 



with ft — \J A 2 + e 2 . The standard expressions for the 
eigenstates of the Hamiltonian dTTI) reduce Eq. ( 10 ) for 
the fidelity to F = 1 - A 2 ^/2f7Vwhere 

rdojS(oj-n) fl-e--/ T e -n/r +e -«/T 



K = 



2ttuj 



T( e -n/T +1 ) 




FIG. 1: The temperature-dependent factor k in the expression 
(JTrj) for the ground state fidelity of an individual qubit in the 
presence of Ohmic environment with cut-off frequency uj c . 



(10) 

We see that the same noise spectral density that defines 



the relaxation and dephasing rates ( 13 ) and ( 14 1 of the 



qubits in the GMQC determines the reduction of the 
ground-state fidelity in AQC. In this respect, the main 
difference between the reduction of fidelity and the real- 
time relaxation and dephasing is that even in the lowest- 
order perturbation theory, the fidelity is reduced by the 
whole spectrum of environmental excitations, and not 
just by limited spectral groups resonant with the qubit 
energy differences or the low-frequency excitations, as in 
Eqs. (pi) and (pi. 



To strengthen this comparison, we consider an Ohmic 
bath characterized by the noise 



S{w) = W(l -e~" /T ) 



(16) 



with cutoff frequency lo c . In this case, the relaxation 



time is T, 



(11) expressed as 



r)(A 2 /n)coth{Q/2T) and the fidelity is 



F = 1 



k 

Q 



k , n 

— tanh — — 
2ry 2T 



(17) 



where Q = T\Q, is the qubit quality factor due to relax- 



ation. Equation ( 15 ) gives the following expressions for 



the factor k at low and high temperatures: 

1 + TT 2 T 2 /3fl 2 , T < Q, 



l_ fln(u; c /fi) 
K - 4tt \0VT) 2 ln( Wc /T), 



t» n. 



(18) 



Equation (17) relates the ground state fidelity to the 
qubit quality factor, Q, as calculated due to relaxation 
only. This shows that the fidelity can be related more 
closely to the relaxation (Ti) and not dephasing (T v ) 
processes. Adding a 1// low-frequency noise of a realistic 
magnitude does not change this conclusion (as discussed 
in more details in the numerical examples below) . As ex- 
pected, a larger Q leads to a better ground state fidelity. 
Figure [l] shows the factor k in Eq. (17) as a function of 
temperature for different cut-off frequencies uj c . It ex- 
(15) hibits the non-monotonic T-dependence, and only weak, 



logarithmic, dependence on uj c , which allows one to es- 
timate the fidelity without precisely specifying uj c . The 
factor k is maximal around fc max ~ 0.5 at T ~ 0.511, 
which leads to a minimum fidelity F ~ 1 — (2Q) _1 . No- 
tice that even a qubit quality factor as low as Q = 10, 
which is practically useless for GMQC, leads to F > 95% 
ground state fidelity. 

We now consider multi- qubit systems, starting with a 
system of TV uncoupled qubits. In this case, the trace in 
the definition of fidelity ^ can be taken independently 
over separate qubits, so that the total fidelity F is the 
product of fidelities Fj, j = 1, ...,N of the individual 
qubits: F = YljFj. For instance, a typical starting 
point of AQC algorithms is to initialize the system in 
the ground state of the Hamiltonian Hjj ( 20 1 . Then, the 



state of all qubits is the same and can be characterized 
by the same fidelity fl7|). Then, 



F = (1 - k/Q) 



N 



-kN/Q 



Q»fc 



(19) 



For independent qubits, fidelity scales exponentially with 
N as a result of the exponential scaling of the probability 
for all qubits to remain in their corresponding ground 
states. Since Q is inversely proportional to the noise 
strength r), by decreasing the noise by a factors of, e.g., 
10, one can achieve the same fidelity with 10 times more 
qubits. 

Next, we focus on how the ground state fidelity behaves 
in practical AQC systems. We use as an example the 
D- Wave One quantum annealing processor installed at 
the University of Southern California (see Ref . [BJ . The 
Hamiltonian implemented by the processor has the form 
of Eq. 0, with 



JY 



N 



N 

£ 



J4. 



(20) 



where hi and Jij are tunable dimensionless bias and cou- 
pling coefficients. The parameters A(s) and B(s) for this 
processor are plotted in Fig. [2)3. We calculate fidelity 
of the ground state for a ferromagnetic chain (illustrated 
in Fig. [2^) with hi=0 and Ji t i+i= — 1, otherwise known 
as a quantum Ising model in a transverse field. Here, 
the length of the chain is varied from A^=2 to 16. Al- 
though this model is exactly solvable (see, e.g., Ref. I2TJ1 
and references therein), fidelity cannot be calculated ex- 
actly for practical noise models in which the coupling to 
environment is dominated by the <t| terms. Hence, we 
calculate the fidelity numerically. In the limit TV— >oo, 
the model is known to have a quantum critical point 
at A(s)= B(s). At this point, the chain goes through 
a quantum phase transition between paramagnetic and 
ferromagnetic phases. In the ferromagnetic phase, the 
ground state is doubly degenerate with respect to simul- 
taneous change of signs of all of terms. Figure [2J: plots 
several of the lowest energy levels of a 10-qubit chain 
relative to the ground state energy E . As seen in this 
plot, the quantum critical point manifests itself as the 
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FIG. 2: a. A ferromagnetic spin chain with transverse field 
and coupling energies given, respectively, by A(s) and B(s) 
in Eq. |TJ. b. Energy scales A(s) and B(s) extracted from 
experimental parameters, c. The lowest 20 energy levels, 
relative to the ground state, of a 10-qubit ferromagnetic chain 
with Jij = —1, as a function of the normalized time s. d. 
Ground state fidelity of the 10-qubit chain of c at T = 20 
mK. The vertical (red) dashed line marks the quantum critical 
point. 



appearance of the doubly-degenerate ground state and 
the minimum in the energy gap between the ground and 
the second excited states. 

To calculate fidelity for this system, we use a realis- 
tic noise model relevant to the D-Wave qubits [IS]. In 
this case, the dominant environmental coupling is to the 
magnetic flux noise, which couples directly to the qubit 
computational basis states represented by the <r| opera- 
tors. The noise spectral density S(u>) was characterized 
in the earlier experiments, which were consistent with the 
noise being a combination of the 1// low-frequency noise 
and an Ohmic noise at high frequencies [19] . For calcula- 
tions of fidelity, we take S(tu) — k(s)[Shf( uj ) + Slp(uj)], 
where Shf(w) is the Ohmic spectral density (16) and 
S lf {lo) = 7 2 /|w|. The coefficient k(s) = B(s)jB(s m ) 
appears because the strength of coupling to flux noise de- 
pends on the persistent current of the flux qubits which 
changes as a function of s (see SI). Here s m is the bias 
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point at which the measurements of r\ and 7 are per- 
formed. Based on the experimental data, we use 77 = 0.1, 
7 = 20 MHz and s m = 0.636. We also assume uj c = 100 
GHz for the high-frequency cutoff and ojl = 1 MHz for 
the low-frequency cutoff (based on a i/~l fis evolution 
time of an algorithm). We found that for these param- 
eters, the fidelity is dominantly determined by the high- 
frequency Ohmic noise and not by the 1/ / noise. 

In principle, since the total number of energy levels 
grows exponentially with N, the time required for nu- 
merical calculation of F also grows exponentially. For- 
tunately, the value of F converges rapidly for a finite 
number of retained energy states. Here, we keep all en- 
ergy levels for N < 10, and up to 2000 energy levels for 
larger chains. The fidelity of the 10-qubit chain is plotted 
as a function of s in Fig. [jjji. The ground state fidelities 
of chains with other lengths (and coupled systems other 
than chains) are qualitatively the same as the one plotted 
in Fig. [2jl. It is clear from the figure that the fidelity is 
minimum close to the critical point. Notice also that the 
fidelity approaches F = 1 as s increases, which is the re- 
sult of Hp commuting with Hj, with only <t| terms and 
negligible other types of coupling to environment. This 
again reflects the fact that the fidelity depends rather on 
relaxation than dephasing. 

Figure [3] shows the numerical results for the ground 
state fidelity for iV-qubit chains with N = 1 to 16 at the 
critical point. For all chain lengths, the fidelity is better 
than 90%. It should be emphasized that these are the 
minimum fidelities at the quantum critical point. The 
fidelity at all other points is larger, and near s = 1, is 
very close to 1 as shown in Fig. [2ji. We have also plot- 
ted in Fig. [3] the ground state fidelity of N uncoupled 
qubits at different N based on the exponential scaling of 
Eq. (19). The scaling and magnitude of the fidelity at 



large N is better for the ferromagnetic chain than for the 
uncoupled qubits. Unfortunately, it was not possible to 
pursue numerical calculations beyond 16 qubits, as direct 
perturbation approximation would break down when F 
strongly deviates from unity. A naive exponential ex- 
trapolation of the data points to N = 128 (representing 
the worse case) still yields F = 0.47, meaning that the 
eigenstates could retain their quantum properties with- 
out error correction for such a large-size system. As in 
uncoupled qubits, if one can reduce the noise by a large 
factor, the size of the chain can be increased by the same 
factor while keeping the fidelity unchanged. In addition, 
other techniques such as dynamical decoupling [21] or 
error correction [33] could be employed to enhance the 
ground state fidelity at large scales. 

Finally, we discuss how the ground state fidelity should 
affect the performance of AQC. In universal AQC Oil], 
the fidelity of the final ground state determines the qual- 
ity of the computation. Indeed, deviations of F from 1 
mean that the statistics of measurements done on this 
state will be different from the one that corresponds to 
the ideal ground state. For instance, in the case of one 
qubit with e = and the Hamiltonian ( 11 ), measurement 
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FIG. 3: Ground state fidelity at the quantum critical point 
for ferromagnetic chains with TV = 1 to 16, at T = 20 mK. 
Circles are the numerical results using (10 1. The red dashed 



line is fidelity of uncoupled qubits from (19 1, with k — 0.32 
and Q = 38.4. 



of a x has a non-vanishing probability 1 — F 2 of produc- 
ing the result a x = — 1 different from the ground state 
(j x = 1 even at temperatures T -C A. However, this ef- 
fect is absent in the special case when the coupling to 
environment via Hj commutes with the final Hamilto- 
nian Hp, leading to F = 1 at the end of evolution, as in 
the adiabatic quantum optimization discussed above and 
shown in Fig. [2ji. In this case, thermal transitions in- 
crease the loss of probability due to small fidelity in the 
middle of the evolution, thereby decreasing the ground 
state probability even further (Pq = PqF 2 ). Therefore, 
the probability will be distributed among the low energy 
states even more than implied by the thermal equilib- 
rium. Part of the probability can be regained later when 
the gap is larger and F is closer to 1. However, since 
the relaxation time becomes exponentially long near the 
end of evolution, the majority of the probability that 
is lost may not be gained back, thus leading to smaller 
probability of success. This makes it important to main- 
tain fidelity close to unity throughout the evolution. We 
stress that most treatments of AQC based on the weak 
coupling master equation, e.g., [5] SOUS], do not take into 
account the effect of deformation of the eigenstates that 
is captured by our calculation of the fidelity. 

In summary, we have proposed using ground state fi- 
delity as a quantity for measuring the strength of deco- 
herence effects in AQC. Fidelity plays a role similar to 
decoherence time in GMQC, but takes into account qual- 
itatively different effects of environment on the ground 
state relevant to AQC. The fidelity is related to the re- 
laxation processes and is relatively insensitive to the de- 
phasing. Our numerical calculations indicate that fidelity 
close to unity can be achieved with a moderate qubit 
quality factor, even for large numbers of qubits. Ground 
state fidelity should be a useful measure of the environ- 
ment related quality of AQC systems in the context of 
further work on important topics in AQC such as quan- 
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turn error correction or the threshold theorem. D.A. Lidar, R. Liu, T. Mahon, F. Nori, A. Yu. Smirnov, 
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SUPPLEMENTARY INFORMATION 

In this document, we provide details of the derivation 
of the main equation for the ground state fidelity and a 
brief description of the numerical calculations. 

I. GROUND STATE FIDELITY 

To be consistent with the notation in the main text, 
we use symbols with (without) "~" to denote quan- 
tities related to the coupled (uncoupled) qubit system 
and environment. We also use m, n, k to enumerate 
qubit system's eigenstates and eigenvalues (e.g., \n), 
E n ). Similarly, the letters v, p correspond to environ- 
ment eigenstates and eigenvalues, and letters a, b corre- 
spond to total system (qubits+environment) eigenstates 
and eigenvalues. The total Hamiltonian is written as 
H = Hs + Hb + Hi, which includes the qubit system 
Hamiltonian 

H s = Y J E n\n){n\, (21) 

n 

and interaction Hamiltonian 

where a" are the Pauli matrices, a — x,y, z, and are 
the noise operators dependent on the heat bath variables 
coupled to the j-th qubit. The bath Hamiltonian Hb is 
not specified, but the bath is assumed to be in thermal 
equilibrium at temperature T. For the purpose of pertur- 
bative calculations of fidelity in this work, the properties 
of q" are completely characterized by its spectral density 

/oo 
dte^iqf^qfiO)) (23) 
-oo 

where qf(t) — e tH Bt e~ lH Bt , and vanishing average 

<#(*)> = o. 

Even if the dissipative environment does not excite the 
system out of its ground state, the non- vanishing Hi term 
changes the structure of the ground state wavefunction 
in comparison with the unperturbed ground state of the 
Hamiltonian Hs ■ The magnitude of this change can be 
characterized quantitatively by fidelity. For two general 
density matrices p and p' acting on the same Hilbcrt 
space, the fidelity (sometimes called Uhlmann's fidelity) 
is given by 

F(p,p')=Try/ y /p ( / y /p. (24) 

If one of the states is a pure state, e.g., p — \ip){ip\, then 
we can use the property ^fp = p = p 2 to write 

= m P p'}} 1/2 - (25) 



Let | a) and |o) denote the unperturbed and perturbed 
states of the total system, respectively. The total density 
matrix of the combined system is p SB = -P a |a) (a| , 

where P a is the probability of finding the total system in 
the state \a). We define the reduced density matrix as 
ps = TtbIpsb], where the partial trace is taken only over 
the environmental degrees of freedom. We also define the 
pure state density matrix p a — |0)(0|, where |0) is the 
ground state of the (isolated) qubit system Hamiltonian 
Hs- Then, we express the fidelity T between p = |0)(0| 
and ps as 

T(p ,ps) 2 = (0\ps\0) - Po , (26) 

where Po is the probability of finding the open system in 
state |0). In equilibrium, this probability is determined 
by the thermal distribution over the system eigenstates 
and also by the effects of the non-vanishing coupling to 
the environment. To separate these two contributions, 
we define the ground state fidelity as 

F^p- 1/2 F{p^ps) = \ (27) 

V M) 

where Po = e~ E °/ T /Z is the Boltzmann probability 
of the ground state |0) with partition function Z = 
e~ £ "/ T in the case of an isolated qubit system. 



II. PERTURB ATIVE CALCULATION OF 
FIDELITY 

In the situation relevant to quantum computation, the 
coupling Hi is weak and can be treated by perturbation 
theory around the non-interacting state of the qubit sys- 
tem and environment. In the context of AQC, one can 
also assume that both the qubits and environment are 
in equilibrium at the same temperature T. These are 
characterized by the density matrices ps = ^2 n Pn\n)(n\ 
and pb = ^2, v Pv\ v )( v \ without the interaction, while the 
total interacting system has the density matrix p S B — 
^ a P |a)(a|, where P a , P„, and P„ are the Boltzmann 
probabilities. The probability P n of finding the qubit 
system in the state \n) in the presence of interaction can 
be written as 

Pn = Tr Bt s[\n)(n\p S B}. (28) 

The interaction Hi creates corrections to the equilibrium 
probabilities P a = P a + SP a , where P a = P„P„, and 
to the wavefunctions: \a(n,v)) = \a(n, v)) + \Sa(n, v)), 
where \a(n,v)) = \n)®\v). In the lowest non-vanishing 
order in Hi, one can separate these corrections in the 
total density matrix of the system: 

Psb = J2i SP -\ a )( a \ + p a\a)(a\] ■ (29) 

a 
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This expression reduces Eq. ( 28 1 for the probability P n 
to: 

P n = SP n +^P m Tr BtS [\n)(n\^p B -\^m,v))(a(m,v)\] . 

rn 

(30) 

Using the relation |n)(n| = 1 — Efc^o *° transform 

the n = in term in Eq. (|30|), we obtain 



Pn — Pn ~\~ &Pn ^ (J- mnPn ^ mnPm) i (31) 
m^n 



Now, we relate the changes in the occupation probabil- 
ities of the qubit states to the energy changes of the total 
system. Expanding the Boltzmann probabilities, we get 

SP a = ~/ C 0E b +P ^ 



-pP a 5E a + f3P a PbSE b 

b 

/3P a Y,(Pb~S ab )SE b . 



(37) 



where 



(n\Tr B [\6a(m, v ))(5a(m, v)\p B ]\n) 



(32) 



The T- factors in ( 32 ) represent the fractions of the prob- 



ability transfer from the state \m) to the state \n) of the 
qubit system due renormalization of the qubit system 
wavefunctions by non-vanishing coupling to the environ- 
ment. Equation (31 1 shows that, in addition to these 



transfers, environment-induced change of the probabili- 
ties P n is also caused by the change in the equilibrium 
probability 8P n due to renormalization of the energy 
eigenvalues (Lamb shift). 

Next, we calculate these two contributions explicitly 
using the perturbation expansion. We start with 5P n . 
As shown below, the changes of the equilibrium occu- 
pation probabilities 5P n of the qubits states are deter- 
mined by the changes 8E a of the total system energies 
averaged over the equilibrium state of the environment. 
Since (q") = 0, the linear corrections to the qubit sys- 
tem energies vanish, and the relevant energy changes are 
given by second-order perturbation: 



SE a = 



(a\Hj 



r-f E a — E b 



E 



j,nm I 



2 IH<^)I 2 



E,j — E.. 



(33) 

The prime sign over the sum excludes terms with zero 
denominator. The average change of the qubit system 
energy eigenvalues due to coupling to the environment, 
which determines SP n , is therefore 



6E n = J2 P^E niV = P v 



E u — Eu 



We express this result in terms of the noise spectral den- 
sities. In this calculation, we assume that the noises with 
different a and j are uncorrelated. The standard spectral 
decomposition, 

s?M = 2^ p,|Hg>)| 2 5^ + £,-£,,), (35) 



transforms Eq. ( 34 1 into 



SE n 



E 

j,at-,m ' 



du lo- 



ot I 2 QCY 



2ir u) n 



ui 



(36) 



Using P a = P n P u , P b = P m P n , and 5 ab = 8 nm 8 ua , we 
see that SP n is indeed determined by the energy shifts 
averaged over the environment: 

sp n = J2 sp - = pJ2 PnP » E( p »^ - & 

)8E m , u 



/3P n J2(Pm ~ S nm )SE m . 



(38) 



Combining Eqs. (36 1 and (38), we obtain 



5P n = -pP n £ (Pm-5 mn )K mk \ 2 [ 



du Sf(w) 

2lT LO km + UJ 



(39) 

Next, we calculate the probability transfer fractions 
r„ m ( 32 ) . It is convenient to view the relation between 



the states with and without interaction as arising from 
the adiabatic switching on of the interaction: 



|2) = U\a) , 



where 



f/ = Texpj-i J dte xt Hj(t)^ 



(40) 



(41) 



Here A is a small positive number which ensures an adi- 
abatic increase of the coupling, and 



E < 



am e^qf(t)\n)(m\ (42) 



with Hj given by Eq. ( 22 ) 



Using Eq. (40) together with (41) limited to the lowest 



order in Hj in Eq. ( 32 ) , we obtain 

dt / dt'Tr B [(n\e M : -H>(i)|m) p B (m\e Ar H^t'^n)) 



mn = Tr B [(n\U\m)p B (m\U^\n)] 



oo 



f dt f dt'Y^Knmfe 

00 J —°B j a 



^ t+t '^^ t - t '\qf(t / )qf(t)}. 
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We combine this expression with Eq. (23) for spectral 
densities to get 



Tmn — ^ Wtnm\ 2 J 2^S?(u) 



3,ot 

-oo J — oo 

S?(w) 



Vl/ |2 /" 



A 2 + (cj nm + cj) 2 



(43) 



In the limit A — > 0, this gives 
The final result is 



duj S?(w) 



(44) 



Pn=Pn-f3P n Krnkl 2 f ~ 

A ... f- W 



27T (W„ m + w) 2 ' 

du Sf{u)(P m ~5 mn ) 



2TT 



(0J mn +U>) 2 



(45) 



This equation gives the normalized fidelity of the 
ground state n — as 



E h%l 2 / 



dw5^(w)(P„-5 n o) 



47T W r , 



■ w 



j,a,n>0 



du;5j'H-(i J n/.Po)S J ? (-^) 

47T (0J n0 + LO) 2 



(46) 



This is our central equation, which is used in the numer- 
ical calculations. As shown in the next subsection, it can 
also be obtained directly from the partition function of 
the total system. However, the physical interpretation 
of the following derivation is less transparent than the 
derivation presented above. 



A. Alternative derivation using partition function 

Consider the partition function 

Z : Tv ( e-^ ) . (47) 



Using Eq. (21 1 one can express the occupation probabil- 



ities P n directly through Z 



1 d 



(3dE T 



InZ = ^Tr 



5 dH 



= Tr (ps B \n)(n\) 



Z \ dE n / 
(n\Tr B (p SB )\n) = (n\p s \n) = P n (48) 



where 



Ps B 



Z 



(49) 



is the total density matrix of the system plus environ- 
ment. 

First, we calculate the partition function in the inter- 
action representation, 



Texp / drHiir) 



(50) 



where H = H$ + H B is the non-interacting part of the 
Hamiltonian and Hj(t) = e^ H Hje~^ H . Therefore, 



Z = Tr e 



= Tr 



- pH Texp / drH^T 



ZTv 



Ps B fTexp J drJJj(r) 



Z( Texp / (1tHi(t 



(51) 



where (...) = Tr[ps B ---] and ps B — e~P H jZ is the den- 
sity matrix for the non-interacting qubit system plus 
environment, which has the form pg B — ps ® p B = 
X^n Pn\n){n\ (g> p B , where P n is the occupation proba- 
bility of state \n) and p B is the density matrix of the 
environment alone. Expanding to second order, we have 

Z = z(l + dr jf dr' {H^H^r'))^ 
= z(l+ f dr f dr'(qf(T)qf(r')) B 



5^ P n (n\a](r)^(r')\r 



(52) 



where (•■■} B = Tr B [p B ...]. We write the integral in this 
expression in terms of the environment spectral density 
using 

{ q «{T) q «{T>)) B = J gs»e-^- '). (53) 
Defining af nm = (n|<r||m) and u) nm = E n - E m , we find 
Integral = £ / f^Knml'W 



dr / dT'e^™-^-^) 



E 



%-P n \a? tnrn \*S?{u>) 



= E 
- E 



(w nm -w) 2 
2?r (w nm -cj) 2 J J 



27T Wnm—W 



(54) 
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In the last step we used P m = e^ Unm P n , which is valid for 
the Boltzmann distribution, and 5?(— £*;) = 5"(w)e _ ^". 
The first integral in the last equation vanishes, giving 



We use M„ m = M mn and change w to — w in some inte- 
grals to get 



Z = Z\ ! 

Therefore, 
In Z = In Z 

and 

P, = 



E 



E 



2-7T w nm -w 



(55) 



F 2 = 1 - p M ™« / 



dw S(ui)P n 



m,n>0 



2-7T W„ 



■ W 



n>0 J 



dw P S(w) - P n S(-u) 



dujPPnK nm \ 2 Sf(u;) 



2tt 



-MP,-1)M [ 



oo 



w n0 + w 
dw S(w 

27T w 



E M «° 



n>0 



1 9 

d 
dE~ 



hxZ 



E 



1 d 



InZ 



2tt 



(56) 



Finally, we transform this equation using the relations 
1 9 



InZ 



PdE n 
d 1 

S-Bn Wfc m +W 
8P m 



Pn 



dE n 



to obtain 



P n =Pn-P J2 K 



mk I 



(w fem +w) 2 ' 

pP„(P m -S mn ), (57) 
dw S"(u>)P n (P m —S mn ) 



E K« 



2tt 

dwP„S?(u/)-P m S?(-u/) 



2tt 



(w mn +w) 2 



(58) 



This is the same equation as (45). Using (27), one recov- 
ers Eq. ( 46 ) for fidelity. 



III. NUMERICAL CALCULATIONS 

In this section, we derive the equations that form the 
basis of our numerical calculation of the ground state 
fidelity F. We assume that the noise is coupled only to 
tr| operators and Sj(u)) — S(ui) is the same for all qubits. 
We also define M mn = lo- 
be rewritten as 



Equation (46) can 



F 2 = l-(3 M mn 
m,n>0 



dw S(u))P n 

2"7T W mn + W 



pYm 0u [^ s( " )Pn 

f- ' J 2tT W „ + W 



n>0 



-Pj2 M ™o 



- M n0 
n>0 



duj S(w)(P -l) 
27T w m0 + w 

do; S(w) - (P n /P Q )S(-w) 
2n (w„ + w) 2 



(59) 



dw S(w)[l - /3(w„ + w)] - (P„/P )S(-w) 



2tt 



(w n0 + 



(60) 



Next, we symmetrize this equation by using the fact 
that£ m>0 P m = l-P : 



P 2 = W E 

™ -v. ^ n « 



do; S(uj)P n — S(—io)P n 



m>n>0 



■| £ P m (M„ 



M 00 ) 



m>0 



<S(w) — 5 f (— w) 

27T W 



m>0 17 

E^/g 

m>0 17 



dw S(w)[l - /3(w m0 +w)] - (P ro /P )5(-w) 



(w m0 + w) 2 
dw /3(w m0 + w)[P 5(w) - P m S(-u)] 



(w m0 + w) ; 



(61) 



Note that each line in Eq. (61) has the form S'(w)/w 



at large w and hence diverges together with the high- 
frequency cutoff in the used model of environmental noise 
(see below). When all the states of the qubit system are 
included, these divergences cancel each other. However, 
only part of the eigenstates are used in the numerical 
calculations, so the divergence does not vanish if one em- 
ploys directly Eq. <[6lj). 



Using the property J2 n l(l n l cr |l m )| 2 = h ^ can be 
shown that 

P m (M mm - Moo) = J2 Mm °(l - Po - Pm) 



m>0 



m>0 



- J2 ( P m+Pn)M„ 



(62) 



m>n>0 



Combining (61) and (62), we get the final result for 
fidelity: 

or, in our perturbation approximation, 

1 



F = I - ~ / dw/(w) , 



(63) 



with 



/(w) =^ ^ M mn f mn {u) + ^ M mQ g m (u 



m>n>0 



m>0 
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Here 



where 



S(h;)P n - S(-o;)P m P m + P n 

w mn + LU 2 

S(u) - S(-u) 
10 

[1 + {P -l)(i{iO ma +Lo))S{uj) 

(uj m0 +uj) 2 
[1 + P p(u m0 +Lj)](P m /P )S(-u) 

i-p -p m sm-sh 

2 uj 



One can see that 



f mn {u>) -> P n S{uj) + 



as w ^ —uj mn , and 



P m + P„ S(w) + S(-w) 



3m (w) -> /3 2 (P ~l/2)SH+£ 



1 - P - Pn S(W) - 5(-w) 



as — > — w m o. This implies that the poles are removed 
and the S(u)/ui divergence at infinity also disappears for 
each term, thus improving the convergence properties of 
this expression for the numerical evaluation of fidelity. 



A. Noise spectral density 

For numerical calculations in this paper, we consider 
the flux noise in aa rf- SQUID. The noise operator q = 
I p 5<fr x is therefore related to the flux noise 5& x through 
the rf-SQUID, where I p is the persistent current of the 
rf-SQUID. The noise is characterized by its correlation 
function through the spectral density 



S(u>) = 



dt 



(q(t)q(0)) = 



dt e iut (S^ lx (t)S^ lx (0)) (66) 



is the spectral density of the flux noise. Since the actual 
noise is of the form of flux noise, the spectral density, 
S(lu) oc I 2 , should depend on the bias point through I p . 
Then, 



S(u) 



M 



r/uje 



\ _ g— /iw/fesT 



with 



(67) 



(68) 



where I pm is the maximum value for I„, r\ is a dimen- 
sionless coefficient characterizing the Ohmic noise and 7 
is an energy scale characterizing the 1/ / noise. 

In the actual quantum annealing system described by 
Eq. (1) of the main text, the persistent current I p and the 
tunneling amplitude of qubits are time-dependent. As 
described in Ref. 18 of the main text, the energy scale 
B(s) depends on I p as B(s) — M I 2 , where M is the 
maximum mutual inductance between the qubits through 
the tunable couplers and s = t/tf. Therefore, one can 
write k(s) = B(s)/B(s m ), where B(s m ) = M I pm - It is 
therefore sufficient to determine s m to calculate n(s) at 
all s. 

Measurement of the 1// noise reveals a flux noise spec- 
tral density of the form S , $(2tt/) = A 2 /\f\ a , with awl. 
From ( 65 1 we find 



7 W V2TTl pm A. 



(69) 



For the qubits in the current D-Wave processors we have 
0.1, A fa 3/i$o, an d I pm ~ 1 MA at s m — 0.636, 



(65) which gives 7//1 ps 23 MHz « 0.02 GHz. 



